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Abstract. The relation between the cosmological evolu- 
tion of the jet-disk symbiosis model for GRBs and the 
cosmic rays energy distribution is presented. We used two 
different Star Formation Rates (SFR) as a function of red- 
shift and a Luminosity Function (LF) distribution to ob- 
tain the distribution in fluencc of GRBs in our model and 
compare it with the data. We show a good agreement be- 
tween the fluence distribution we obtain and the corrected 
data for the 4B BATSE catalogue. The results we obtain 
are generally valid for models that use jet physics to ex- 
plain GRB properties. 

The fluence in the gamma ray band has been used 
to calculate the energy in cosmic rays both in our Galaxy 
and at extragalactic distances as a function of the redshift. 
This energy input has been compared with the Galactic 
and extragalactic spectrum of cosmic rays and neutrinos. 
Using our jet disk symbiosis model, we found that in both 
cases GRBs cannot give any significant contribution to 
cosmic rays. We also estimate the neutrino background, 
obtaining a very low predicted flux. We also show that the 
fit of our model with the corrected fluence distribution of 
GRBs gives strong constraints of the star formation rate 
as a function of the redshift. 

Key words: Gamma ray : bursts - Stars : statistics - 
ISM : cosmic rays - Cosmology : distance scale 



1. Introduction 

More than 30 years after their discovery, thanks to 
the Burst and Transient Source Experiment (BATSE) 
and the Italian-Dutch satellite BeppoSAX, the scien- 
tific community now knows that Gamma Ray Bursts 
(GRBs) are isotropically distributed in the sky (Fish- 
man & Meegan 1995) and that at least some of them 
are at cosmological distances (GRB970228: Djorgovski et 
al. 1999b, GRB970508: Metzger et al. 1997, GRB971214: 
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Kulkarni et al. 1998, GRB980613: Djorgovski et al. 
1999a, GRB980703: Djorgovski et al. 1998, GRB990123: 
Hjorth et al. 1999, GRB990510: Vreeswijk et al. 1999, 
GRB990712: Galama et al. 1999). But the present data 
available for redshift position and host galaxy localiza- 
tion are still too few to give us good statistics to study 
the evolution of GRBs and their redshift distribution. Be- 
fore the discovery of GRB afterglows by BeppoSAX, the 
only way to study their distributions was to compare some 
GRB properties (like for example the intensity) , with some 
parametric models (Fenimore & Bloom 1995, Cohen & Pi- 
ran 1995, Kommers et al. 1999). Because of this lack of 
information, it is still necessary to assume that GRBs fol- 
low the statistical distribution of some other better known 
objects to obtain the GRBs fluence or flux distribution it- 
self. 

The origin of GRBs is still controversial. According to 
different models, their progenitor can be identified with 
the merging of two neutron stars, or with the collapse of 
a massive star. In the model presented by Pugliese et al. 
(1999), GRBs are created inside a pre-existing jet in a bi- 
nary system formed by a neutron star and an O/B/WR 
companion, where the input energy comes from the col- 
lapse of the neutron star into a black hole and the emis- 
sion is due to synchrotron radiation from the ultrarela- 
tivistic shock waves that propagate along the jet with a 
low-energy cut-off in the electron distribution. Following 
this scenario, the birth of GRBs cannot happen too far 
from the region where the progenitor formed, and this 
implies that their rate should be connected with the Star 
Formation Rate (SFR). Already other authors studied the 
connection between the SFR and GRBs flux distribution. 
For example, Wijers et al. (1998) showed that the assump- 
tion that the GRB rate is proportional to the SFR in the 
universe is consistent with the GRB flux distribution. 

In Sect. 2 we calculate the cumulative distribution of 
GRB fluences using two SFR distributions as a function 
of redshift, the one by Miyaji (Miyaji et al. 1998), and 
the other by Madau (Madau et al. 1996). We compared it 
with the data from the BATSE catalogue. In Sect. 3 we 
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calculate the maximum energy available in our model to 
obtain high energy cosmic rays. In Sect. 4 we present our 
results for the contribution of GRBs to the cosmic rays 
distribution, both Galactic and extragalactic and in Sect. 
5 the eventual contribution from GRBs to the neutrino 
flux. 

2. The rate of GRBs 

The aim of this section is the calculation of the number of 
GRBs per year and per 100 Mpc 3 as a function of redshift, 
assuming beamed emission in a jet, (e.g. Pugliese et al., 
1999). The following calculations arc quite general and 
we emphasize that their validity does not depend on the 
particular model used. 

To obtain a result as close as possible to the data, we 
chose the fiuence as the quantity that can well represent 
the characteristics of GRBs. We follow Petrosian & Lloyd 
(1997), who showed that the fiuence is the most appro- 
priate parameter to study the cosmological evolution of 
GRBs. The relation between the redshift z and the fiu- 
ence / is given by (Petrosian & Lloyd 1997) as a function 
of the specific luminosity L expressed in [erg s _1 Hz -1 ]: 



/ = 



L At Av 



erg/cm 2 . 



(1) 



Here a is the photon flux spectral index, = (2c/h)(l + 
z — ^/l + z) is the cosmological luminosity distance (Wein- 
berg, 1972), and H Q = h (100 km s -1 Mpc -1 ) is the Hub- 
ble constant. In the transformation from the emitter frame 
to the observer frame the two corresponding redshift con- 
tributions from the frequency and the time dependence 
are cancelled. 

In our model (Pugliese et al. 1999), we used a = 2, 
and the corresponding value for the fluence is: 



/ 



LAtAv 



(1 + z) erg/cm 2 . 



(2) 



In this way the redshift as a function of the fluence is: 

l + </,W=Q^+l) , (3) 

where /+ = L ^ 1 4 ^ tA ' y ^r erg/cm 2 is the reference fluence 
and L(vi) is expressed in [erg s _1 Hz -1 ]. 

At this point of our calculations it is important to de- 
fine the role of the parameter In fact a relevant ques- 
tion is whether GRBs are standard candles (i.e. /* is con- 
stant) or whether they are distributed with a Luminosity 
Function (LF) with, e.g., a power law in /*. Developing 
our calculations, we arrived at the same results found by 
Kommers et al. (1999). Here the authors showed that the 
best model is the one in which the star formation rate 
is combined with the luminosity function distribution. In 
agreement with them, also in our model GRBs cannot be 



considered as standard candles. We use a distribution for 
the luminosity function with a power law index — /?, and 
the corresponding law for the parameter /* is given by: 

4>{f*) d f* = <hf* P 4f*. (4) 

where </>o is the normalization parameter equal to l/(/~ b 1 — 
We adopt /* a - 10- 8 erg/cm 2 , and /* b - 2.2 x 
10~ 4 crg/cm 2 . 

Here we do not yet answer the question what the 
physics of this luminosity function may be. It is plausible 
in the context of our model, that it is directly connected 
to the mass flow in the pre-existing jet prior to the GRB 
explosion. If this were the correct interpretation, the mass 
accretion rate, and correspondingly, the mass flow rate in 
the jet may follow a power law distribution, a point which 
we will pursue elsewhere. 

We also calculate the GRB rate using two different 
star formation rates and compare the corresponding re- 
sults with the data. 

2.1. The number count 

The number count of GRBs sources is given by the expres- 
sion dN(z) — F{z) (dt/dz) dV dz, where F(z) is equal to 
the product of the SFR ip{z) and the luminosity function 
4>{z). Obviously the luminosity function may change with 
redshift z, but for simplicity and as a first step, we use this 
ansatz. In term of the fluence / of GRBs in the gamma ray 
band, and of the reference fluence the number count 
is: 



dN(f,u) = foH/,/*)] 

x[dz (/,./*)], 

where: 



dz 



(5) 



• The first term ip(z(f, /*)) represents the star formation 

rate as a function of the redshift. 

• The second term represents the temporal interval in 

which the rate is calculated. It is given by: 

(6) 



, (/>/*) — . ■ 



dz fc (WE77+i) 

The third term is the interval of volume as a function 
of the fluence. It corresponds to the value: 



d 



■dz 



Kl4 = ^7 517 + 1 17 + 1 (7) 



h 3 / V 2 V / 



• The fourth term represents the redshift interval and it 
is given by: 



dz{f,U) 




(8) 



27T (v7) < 

The integration over the parameter /* will be shown 
in the Eq. 11. 
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2.2. The GRB rate distribution using the Miyaji SFR 



Here we use the star formation rate as a function of the 
redshift presented by Miyaji et al. (1998), and compare the 
corresponding GRB rate with the data corrected for selec- 
tion and calibration effects by Petrosian (priv. comm.) for 
the 4B BATSE catalogue. 

In their paper Miyaji et al. (1998), calculated the dis- 
tribution of Seyfert galaxies as a function of redshift. We 
considered this same distribution for the SFR, and ap- 
proximated it with the following function: 



ip(z) = A exp(a z), 



(9) 



where z is the redshift, A ~ 10~ 6 [/i 3 Mpc -3 ], a = 2.5, 
and we used f2 = l,A = 0asa simple reference. This law 
is valid up to z = 2, and then continues as a constant to 
high rcdshifts. 



q_ 3.5 

"O 

CI) ? 
> J 

L2.5 

O 



1.5 



0.5 





















\ '-■ N '. 


" , , , , 1 , , , , 1 , , , , 1 , , , 


\ - \ 



-7.5 



-7 



-6.5 



-5.5 



-5 -4.5 -4 

Fluence log(f) 

Fig. 1. The corrected data (Petrosian priv. comm.) from 
the 4B BATSE catalogue (full circles) are compared with 
the distribution of GRB fluence calculated using the 
Miyaji SFR and a power law luminosity function index 
f3 = 1.5 (dotted line), (3 = 2.0 (dashed line), and (3 = 3.0 
(solid line). 

This SFR as a function of the fluence is given by: 



</>(/,/*) = A exp 




(10) 



Using the LF given in Eq. [|, the rate of GRB fluence 
is obtained from the product of the terms of the Eqs. || 
- 1[ and Eq. where we now integrate over that is: 



N(f) df 



2^AMr df) / *- 1/2 



exp 
1 



2VX+1! -I 



y/x + 1 ( Vx + 1 dx 



(11) 



where x = /*//, and x max and x m i n are the limits of 
integration defined as the maximum and the minimum of 
the intersection between the interval [/* a , /*b] relative to 
the LF distribution and the interval in which the ratio 
/*// is defined. 

In Fig. ^ we plotted the cumulative curves correspond- 
ing to different luminosity function power law indexes 
(j9 = 1.5 dotted line, (3 = 2.0 dashed line, and (3 = 3.0 
solid line) in Eq. [H] and compared them with the corrected 
data of the 4B catalogue from BATSE. It is evident that 
even if we change the luminosity function distribution in- 
dex, we cannot obtain any better fit or calculate any GRB 
rate using this SFR distribution. 

2.3. The GRB rate distribution using the Madau SFR 

In this paragraph we calculate the GRB rate based on the 
same procedure as in the last section, using the SFR as a 
function of redshift presented by Madau et al. (1996) and 
the luminosity function distribution of Eq. ^. 
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Fig. 2. Madau SFR curves at different slopes after the 
redshift z* = 1.7. The full circles are the data points from 
Steidel (1999). In this article, Steidel considers the curve 
ai = 1 as one of the many possible curves consistent with 
both the current data on the far-IR background and the 
galaxies detected in the new submillimiter band. 
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We approximated the Madau SFR by the combination 
of two exponential functions. The first part of the curve 
is: 

fci(z) ~ exp(ai z), (12) 

where z is the redshift, and a\ = 2.9 is a constant given 
by the IR source counts. This law is valid up to z* — 1.7, 
after which it is substituted by the function: 

k 2 (z) ~ exp(— a 2 z), (13) 

and a 2 is a constant. 

Because of the uncertainties of the dust extinction at 
high redshift and also the difficulties in the redshift deter- 
mination for the SCUBA sources (Sanders 1999), the star 
formation history beyond the redshift z = 2 is still un- 
clear, therefore we considered different slopes for this part 
of the Madau SFR curve and plotted them in Fig. [2] to- 
gether with the experimental data given by Steidel (1999). 

The cumulative GRB fluence distribution is obtained 
from Eq. [ll], where the SFR function is substituted with 
Eq. [l2| and Eq. [l^. The parameters that we can change 
to fit the data are the redshift z* , the power law index 
a 2 in the Eq. |l3] and the power law index (3 and the up- 
per limit /*b of the interval in which the luminosity func- 
tion distribution is defined. Following one of the possi- 
ble curves that fit the data shown by Steidel (1999), we 
chose z* = 1.7, (3 determines the slope of the cumulative 
distribution curve, a 2 defines curves with different slope, 
and variations in /*b correspond to little changes in the 
part of the curve relative to the strongest GRBs, we chose 
/*b = 2.2 x 10 _4 erg/cm 2 . To reproduce the corrected data 
by Petrosian of the 4B catalogue, we probed all the differ- 
ent cases changing the slopes of the SFR and the power 
law index of the fluence distribution. The only values for 
the power law index of the luminosity distribution that 
successfully fits the data was (3 = 1.55. 

In Fig. |^ we show the cumulative fluence distribution 
corresponding to the different slopes of Fig. || and the 
power law index f3 = 1.55. We compared these curves with 
the data corrected for selection and calibration effects by 
Petrosian for the 4B BATSE catalogue: we can fit these 
data only if a 2 is in the range 0.8 ~ 1.3. In Fig. |^ the data 
are fitted with a 2 = 1.0. 

Considering the total number of GRBs in the BATSE 
catalogue, an observing time of 8 years, a volume scale of 
/i- 3 10 10 8 Mpc 3 , a beaming factor ^ = 200 6^, wne re 
6*_!j = 6'j/(10 _1 rad) is the opening angle of the jet, and 
the factor 22 coming from the integral of the distribution 
in Eq. [H^ calculated with these new SFR and LF, we ob- 
tain the following rate of GRBs: 

lO-^i^ezlj) GRBs (yr)" 1 (100 Mpc -3 ). (14) 

This value obtained considering beaming effects is 
close to the number given by Piran (1999). 
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Fig. 3. Cumulative curves corresponding to the different 
slopes of Fig. |[ 
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Fig. 4. Comparison between the distribution of GRB flu- 
ence using the power law Luminosity function </>(/*) oc 
f* 1 ' 55 ! an( i the Madau SFR (solid line) and the 4B cor- 
rected data from Petrosian (priv. comm.) (full circles). 

The result depends strongly on the power index, while 
it is not influenced by the interval [/* a , /*b] of the LF dis- 
tribution we used. We defined a lower limit for the fluence 
equal to 10 _8 erg/cm 2 , and a ratio /*b//*a equal to 10 4 . A 
change in the upper limit of this interval corresponds to a 
small change in the tail of the fluence distribution curve, 
for x > 10~ 4 . In our model, this interval in the fluence 
corresponds to an interval in the initial energy deposited 
in the jet in the range [10 48 , 10 52 ] ergs. 

3. Maximum energy available 

Another important step is to check if it is possible to 
produce neutrinos and high energy cosmic rays with our 
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model (Pugliese ct al. 1999). Therefore we calculated the 
maximum energy available in the emission region. Irre- 
spective of any acceleration mechanism, the maximum en- 
ergy of charged particles, here protons, is that given by the 
spatial limit in the comoving frame, where the gyromotion 
just fits the space available: 

E^= 1 - %h e(B i Azf h \ (15) 

Here e is the charge of an electron, B is the mag- 
netic field in the shock frame given by B^ si \t) = 
10.24 (^^^/^(e-V^Z^Vj)*"^, and Az is 
the thickness of the emission region. For the calculation 
of the thickness we have two possibilities, and each of them 
corresponds to one of the following cases: 

• For the first case we used the thickness of the emission 

region in the shock frame (sh), given by z/(47 s h): 

= 7sh (e S (sh) ^) - 1-91 x lO^(El( 4 M%] x 
Vo^o 1/2 0Zl^ 2 2 )t-^ eV. (16) 

• In the second case we considered the width of the shell: 

£"i° x b) = 7sh(e 5 (sh) ^j) * 0.93 x 10 23 i4 /2 x 

{^ 1/2 lZ) t~ 1/2 eV. (17) 

At this point it is necessary to check when the thick- 
ness of the emission region is lower than the width, which 
corresponds to fiS < A b *. To do this we calculated 
the time at which these two quantities are equal: 

He&=e&) ~ 5 ' 62 x ^(^^VsV-ij s- (18) 

It means that for about two months iS'max^ < eJ!!?. 
Obviously, only E^ty} is valid, so after 10 seconds the 
upper limit of the maximum energy available in our model 
is about 1.07 x 10 21 eV, while after two months the energy 
to consider is E^_^ . Adiabatic losses will diminish these 
energies for charged particles like protons. 

4. Cosmic ray contribution 

Cosmic rays are ionized nuclei, mainly protons, that ex- 
tend from low energies (few hundred MeV) up to very 
high energies (about 3 x 10 20 eV). Their spectrum is de- 
scribed by a power law (dN/dE) = E~( K+1 \ and shows 
two breaks in the slope. From low energies up to about 
5 x 10 15 eV, known as the knee, the spectrum follows a 
pure power law with k ~ 1.7. The detailed shape of this 
break and the precise position are still unknown. Beyond 



the knee up to about a second break point at 3 x 10 18 
eV, known as the ankle, the pure power law has an index 
k ~ 2. 

It has been proposed (Biermann 1993, Stanev et al. 
1993), that three components contribute to the cosmic ray 
spectrum: a) explosion of supernovae into a homogeneous 
interstellar medium (ISM), accelerating particles up to en- 
ergies of about 10 5 GeV. The spectrum for these particles 
is a power law with an index of -2.75, after considering 
the leakage from our Galaxy, b) Explosion of stars into 
their former stellar wind (like Wolf Rayet stars), produc- 
ing particles with energies up to about 3 x 10 9 GeV. The 
corresponding spectrum switches at the knee from -2.67 to 
-3.07 and this difference in the spectral index derives from 
a diminution of the particle curvature drift energy gain. 
In Bicrmann's cosmic ray model for the Galactic compo- 
nent (Biermann 1997), the energetic protons are produced 
in the shocks of supernova explosions in the interstellar 
medium, while all the heavier elements are produced in 
the shock waves propagating in the stellar wind of the 
progenitor star, c) Production of particles with energies 
up to 10 12 GeV from the hot spots of Fanaroff Riley class 
II radio galaxies. Their spectrum has an index -2 at the 
source, and one needs to take the interaction with the 
cosmological microwave background into account. 

We expect the spectrum of cosmic rays above about 

5 x 10 19 GeV to be strongly attenuated because of the 
interaction of nuclei and protons with the 2.7 K cos- 
mic microwave background, giving rise to the so-called 
Greiscn-Zatsepin-Kuzmin (GZK) cut-off. The extragalac- 
tic sources cannot produce all the total cosmic ray energy 
density observed at Earth, but they could give a contribu- 
tion to the ultra-high energy (UHE) part of the spectrum. 
The origin of the cosmic rays above 3 x 10 18 eV is not 
yet clear, but the common idea is that they are extra- 
galactic and probably connected with the most powerful 
radio galaxies (Biermann & Strittmatter 1987, Bcrczinsky 

6 Grigor'eva 1988, Rachcn & Biermann 1993, Rachen et 
al. 1993). At the moment, our knowledge of the sources of 
the highest energy cosmic rays is limited by the small num- 
ber of events detected by the present experiments. Con- 
sidering that at 10 20 eV the rate of cosmic rays is about 1 
event per km 2 per century, it is clear that to detect them 
it is necessary to have both large aperture detectors and 
a long exposure time. 

In this general context, it is interesting to check the 
eventual Galactic and extragalactic energetic contribution 
given by GRBs to the cosmic ray spectrum in our model. 
In fact GRBs seem to be very powerful explosions and they 
inject a large amount of energy and elementary particles 
into the interstellar medium. 

We used the GRB rate obtained with the SFR from 
Madau and two different energetic approaches to calcu- 
late the contribution from GRBs to the cosmic rays and 
the neutrino spectra. First we considered that each GRB 
gives the same contribution equal to 10% of a fixed initial 
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energy of 10 51 ergs. Secondly we assume that each GRB 
contributes proportional to its own fluence, which we as- 
sume is distributed with a power law, corresponding to a 
range of initial energy of [10 48 , 10 52 ] ergs. 

4-1- Extragalactic contribution 

For GRBs it is important to identify which particles con- 
tribute to the cosmic ray flux. As Rachen & Meszaros 
(1998) showed, during the main burst protons lose most 
of their energy because of adiabatic expansion, while neu- 
trons can be better candidates to obtain ultra high en- 
ergy cosmic rays (UHECR) and neutrinos. In fact neu- 
trons carry about 80% of the proton energy, and because 
they are not coupled to the magnetic field, they can es- 
cape the fireball and through the /3-decay give a cosmic 
ray proton spectrum. We followed this same logic in our 
calculations below. 

a) We assumed that the total energy discharged by 
each GRB in the ISM for hadrons as well as neu- 
trinos is equal to 10% of the initial energy E 5 i = 
£'/(10 51 erg) deposited in the jet, that is 7710 = 
(10%£)/(10 50 erg). It means that in the Hubble time, 
the energetic contribution per unit volume inside the 
whole universe given by all the GRBs is equal to 
lO- ai -°(/i 3 0I^»7io) erg/cm 3 . 

To derive the spectrum of GRBs and to com- 
pare it with the one of cosmic rays out of our 
Galaxy, we need to calculate the normalization fac- 
tor N in N (-Jr-) dE, where N is expressed in 
[GeV" 1 cm" 2 s^sr" 1 ]. 

N is obtained directly from the integration of this 
power law, remembering that the result of this inte- 
gral is equal to the energy per volume produced by all 
the GRBs: 



lO- 21 - ^:^). (19) 



nE2 












Ie 1 







The corresponding value is 



N 



lO- 8 S (h 3 9zlmo) GeV- 1 cm" 2 s" 1 sr" 1 . 



(20) 



b) We assumed that the total energy discharged by each 
GRB in the ISM for hadrons as well as neutrinos is pro- 
portional to its own fluence in the 7 band. To obtain 
the energetic contribution per unit volume inside the 
whole universe given by all the GRBs in the Hubble 
time, we integrated the Luminosity Function distribu- 
tion and had 10~ 23 - 6 (/i 3 6Cij?7io) erg/cm 3 . 
The corresponding normalization factor N in the curve 
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Fig. 5. Comparison between the extragalactic GRB con- 
tribution for the case a) (solid line) and the case b) 
(dashed line) and the all particle cosmic ray spectrum, ex- 
pressed in [GeV" 1 cm" 2 s" 1 sr" 1 ]. The dotted lines show 
that it not possible to have any contribution beyond 10 18 
eV, because of the interactions with the microwave back- 
ground. The stars represent the cosmic ray data from 
Akeno experiment (Nagano et al. 1984), the open squares 
are the Fly's Eye data (Baltrusaitis et al. 1985) and the 
full circles are the AGASA data (Yoshida et al. 1995). 



In Fig. |5| we plot the all particle energy spectrum as 
measured by different ground-based experiments (Bier- 
mann & Wiebel-Sooth 1999), the spectrum in the case 
in which each GRB gives the same contribution to CR's 
(solid line) and the spectrum corresponding to a contri- 
bution from each GRB proportional to its own fluence 
(dashed line). The dotted lines represent unlikely contri- 
butions because beyond 10 18 eV the interactions with the 
cosmological microwave background are relevant and they 
make the curves much flatter, and therefore much lower. 
From this graphic it is clear that in our model GRBs do 
not give any energetic contribution to the extragalactic 
cosmic ray spectrum. 

4-2. Galactic contribution 

To calculate the contribution from GRBs to the Galactic 
cosmic ray spectrum it is necessary to know the GRBs pro- 
duction rate in our Galaxy. We considered that the ratio 
between the in/out GRB rates is equal to the ratio of the 
infrared (IR) luminosity inside and outside our Galaxy. In 
our Galaxy this luminosity is Lir, — 10 10 Lq ~ 10 43 6 erg/s 
at 60/^m. Using the equation (1) of Malkan & Stecker 
(1998), we have an extragalactic IR luminosity Ltr — 
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10 44 f5 /i 3 erg(100Mpc 3 ) 
GRB rate is equal to 10" 



. This means that the Galactic 
6A (h 3 eZ%Vio) GRBs per year. 



a) Following the same procedure used for the extragalac- 
tic case, we calculated first the case in which each GRB 
gives the same contribution to the cosmic ray spec- 
trum. We obtained a diffuse density energy of GRBs 
m our Galaxy equal to lO _16 -°(/i 3 0Z^io) erg/cm 3 . 
In our Galaxy, cosmic rays have a time scale to escape 
with an energy dependence that goes as E~ x / 3 (Bier- 
mann 1995, Biermann et al. 1995). This term has to be 
taken in account to calculate the total spectrum of the 
primary cosmic rays inside our Galaxy. It means that 
it is necessary to multiply the injection spectrum E~ 2 
times the leakage term E~ x l 3 to obtain the final plot 
of the contribution from GRBs to the Galactic cosmic 
ray spectrum: N (J^) 7 ^dE. The normalization fac- 
tor is obtained using the same procedure of the last 
section: 



f 2 



— n(— 

c \E 



EdE 



= 10 



- 16 -°(/ l 3 ^ j r 7l0 ). (22) 



The corresponding value is 



N ~ lO _3 - 9 (/i 3 0I^i O ) GeV" 1 cm" 2 s" 1 sr" 1 . (23) 
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Fig. 6. Comparison between the Galactic GRB contri- 
bution for the case a) (solid line) and the case b) (dashed 
line) and the all particle cosmic ray spectrum, expressed in 
[GeV" 1 cm -2 s _1 sr -1 ]. The dotted lines and the points 
are the same of Fig. [5| 



b) For the second case we assume that each GRB con- 
tribution is proportional to its own fluence in the 7- 
band, and integrating the luminosity function distri- 
bution we had an energy density in our Galaxy equal 
to 10- 18 7 (/i 3 6'I 2 J .r7io) erg/cm 3 . 

The corresponding normalization factor N in the spec- 



trum N 



dE is 



N ~ lO- 6 - 5 (/i 3 0I 2 1? 7io) GeV" 1 



cm 2 s 1 sr . 



(24) 



In Fig. |^ we compared the all particle energy spec- 
trum as measured by different ground-based experiment 
(Biermann & Wiebel-Sooth 1999), with the spectrum from 
GRBs in the case that each of them gives the same con- 
tribution (solid line) and with the one in which the con- 
tribution is proportional to the fluence (dashed line). The 
dotted lines beyond 10 18 eV show where interactions with 
the microwave background may become relevant. Since 
massive star formation is highest in the Galactic central 
region, any contribution is limited to energies for which 
the Larmor radius becomes as large as the Galactic disk, 
and the AGASA data suggest only a small anisotropy at 
the Galactic center. But even if from an energetic point 
of view, GRBs could give a contribution at these high 
energies, this is ruled out considering that the temporal 
interval between two GRBs in our Galaxy is larger than 
the leakage time of cosmic rays at these energies. There- 
fore, in our model also in our Galaxy GRBs cannot give 
any contribution to the cosmic ray spectrum. 

5. Neutrino production 

We also would like to probe the energetic contribution 
of GRBs to the neutrino flux. Both the very high energy 
(VHE) neutrinos, with energies in the range 10 10 -^10 17 eV, 
and the ultra high energy (UHE) neutrinos, with energies 
E > 10 17 eV originate in the interactions of protons with 
photons, through the reactions P7 — > mt + , n + — > fJ. + ^, 
fi + — > e + v^i7 c . The efficiency for the neutrino production 
depends on the fraction of proton energy converted into 
charged pions and on how much energy pions and muons 
keep before decaying. 

Some authors (Waxman & Bahcall 1997), proposed 
that GRBs can be associated with neutrinos produced 
in the photohadronic reaction. Pions come from the in- 
teractions between accelerated protons and gamma rays 
in the fireball, and their decay produces neutrinos and 
anti-neutrinos together with other elementary particles. 
Waxman & Bahcall (1997) gave an upper limit for the 
corresponding neutrino flux. 

Other authors (Rachen & Meszaros 1998), argued that 
the upper limit obtained in this way is optimistic. In fact, 
protons emitted in the earliest burst do not have enough 
energy to leave the expansion region because of adiabatic 
losses. Instead neutrons can easily escape and contribute 
to the neutrino flux. 
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To check what is the neutrino flux in our model, we 
calculated the initial photon density number in the 
interaction from the synchrotron photons in our model. 
We obtained a number of photons/cm 3 equal to: 



L (sh)/f 47r ^ 

V 4 7sh 



100 KeV 



7sh 



l/4 A >5/4 
5j 



6.94 x 10 lb (£?^i ' M_' 
,7?/ 4 ^- 2 .f-7/4 



J 0.3 



(25) 



To obtain the number of hits in the 737 collisions, we 
used a cross section er p7 = 2.70 x 10~ 28 cm 2 , correspond- 
ing to the maximum of the curve describing this reac- 
tion. This implies a number of hits per proton equal to 

0.99 {E-^M%y,r) » i > ' ■ 



"10 



o 



10 



Freius limit at 90% CL- 



IO 



10 



10 



10 10 10 

log(E/GeV) 



Fig. 7. The GRB contribution to the neutrino spectrum 
from our model, expressed in [GeV -1 cm -2 s _1 sr -1 ]. The 
dashed line represents the measured upper limit to the 
neutrino flux from the Frejus experiment (Rhode et al., 
1996). Beyond 10 9 GeV, corrections from the proton and 
neutron interactions with the microwave background can 
be expected. 

To calculate the energy rate of neutrinos in our model 
we assumed that about 20% of the energy of protons 
E p — 10 50 erg goes into neutrinos and that in a first 
approximation there are no multiple hits for the same 
proton in the interaction. Following the same pro- 
cedure used to calculate the extragalactic cosmic ray 
contribution from GRBs, we obtained that the normal- 
ization factor corresponding to the spectrum (E/Eq)~ 2 
in the case in which each GRB gives the same con- 
tribution to the neutrino spectrum is N ~ 1.55 x 
10~ 10 (/i 3 6> 2 y.rjio) GeV^ 1 cm~ 2 s _1 sr^ 1 . This spectrum is 



plotted in Fig. [?] and it shows that in our jet model GRBs 
give a contribution to the cosmological neutrino flux at a 
low level. 



6. Discussion and conclusions 

In the first part of this article we calculated the GRB rate 
and compared the corresponding cumulative distribution 
in fluence with the observational data. There were two 
main points to decide on: a) which luminosity function 
distribution and b) which star formation rate were the 
best to reproduce the data. 

We checked if in our jet model GRBs were standard 
candles. But we did not obtain any good fit, therefore 
we tried a power law for the luminosity distribution, 
0(/) ■ Together with this function we used the SFR 

from Miyaji et al. (1998), in which the SFR grows linearly 
from z = up to about z = 2 and after this redshift is flat 
up to z = 6. There was not a good agreement between the 
distribution in fluence we obtained and the corrected data 
by Petrosian for the 4B BATSE catalogue. Therefore we 
used the SFR model from Madau et al. (1996), in which 
the rate follows the same behavior of Miyaji's up to a red- 
shift z*, and after this redshift the SFR begins to decrease. 
There are still some uncertainties about the shape of this 
second part of the distribution, so we considered different 
curves with different slopes (depending on a constant 02) 
for this decreasing part. We have only two free param- 
eters that we can change to fit the data, the power law 
index /3, and the exponential index ai. The redshift z*, 
and the upper limit /*b of the interval in which the LF is 
defined are not really free parameter because their range 
is limited by the observations. We can reproduce the 4B 
BATSE corrected data by Petrosian using the following 
values: redshift z* = 1.7, /3 = 1.55, /* b = 2.2 x 10~ 4 and 
ai in the range [0.8,1.3]. These parameters are remark- 
ably constrained. Thus, given a final model for GRBs and 
a cosmological model, we may be able to derive strong 
limits on cosmological parameters. 

The key point of the second part of the article is the 
calculation of the GRB rate inside and outside our Galaxy. 
The results we obtained depends mainly on three other pa- 
rameters, the value of the Hubble constant H , the open- 
ing angle of the jet #j, and the power law index a we 
assumed for the electron distribution in our model (see 
Eq. 0). 

The choice of a = 2 has been done in the first version 
of our model, where even if we simplified in many places 
the physics used, we obtained a good agreement with the 
data. Any possible small changes in this interval will influ- 
ence the results obtained in this work in a marginal way, 
because the dependence on a in the expression of N is not 
strong. However, putting a = 3, which is not suggested by 
the data, would strongly influence our results. 

On the other hand, changes in H Q and 9 3 , because of 
the strong dependence in N will influence the results. It is 
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interesting to note that both, a lower value of the Hubble 
constant and a smaller opening angle of the jet, go into 
the direction of decreasing TV. But these changes cannot 
influence our results because the interactions of cosmic 
rays with the microwave background and the large differ- 
ence between the Galactic GRBs and the leakage time of 
CR's rule out that in our jet model GRBs can give any 
contribution beyond 10 18 eV. 

In the calculation of the contribution from GRBs to 
the neutrino flux, it is important to define the parameters 
that characterize the neutrino production. The number 
of hits in the pry collisions has been obtained considering 
only the major photohadronic interaction channel, the one 
that gives a single pion. At higher energies it is possible 
to have the channels for the production of multi-pions, 
pj — > n27r + 7r~ and p7 — > 7j37r + 27r - , but as the energy in- 
creases the corresponding cross sections decrease. We used 
only the first channel with the highest cross section, be- 
cause the energies involved are not high enough to require 
secondary channels. 

In the context of our jet-disk symbiosis model for 
GRBs, there is only one way to make the extragalactic 
contribution to CRs significant at the highest energies, 
and that is to drastically increase the CR energies per 
GRB deposited. We consider this implausible in the con- 
text of our model because of energy conservation require- 
ments. We use 10% of the entire energy available, and so 
the CR contribution may be increased over our simple cal- 
culation by a factor of a few, but not more. This implies 
that the contributions given by each GRB to CR and neu- 
trino flux cannot be much bigger than the energy emitted 
in the 7-ray band. 

Using a relatively small set of parameters, the jet-disk 
symbiosis model applied to GRBs, a tested star formation 
rate and the fundamental physics of the photohadronic 
interactions we arrive at the conclusion that GRBs are 
not able to give any significant contribution to the high 
energy cosmic ray spectrum both inside and outside our 
Galaxy and predict only a low flux of neutrinos. 

A main conclusion of this work is that fitting the cor- 
rected fluence distribution of GRBs with the jet-model is 
well possible. The fit allows strong constraints of the star 
formation rate as a function of rcdshift. 
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